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We develop a supersymmetric field theoretical description of the Gaussian en- 
semble of the almost diagonal Hermitian Random Matrices. The matrices have 
independent random entries Hi>j with parametrically small off-diagonal elements 
Hij/Ha ~B< 1. We derive a regular virial expansion of correlation functions in 
the number of "interacting" supermatrices associated with different sites in the real 
space and demonstrate that the perturbation theory constructed in this way is con- 
trolled by a small parameter B. General form of the integral expression for the m-th 
virial coefficient governed by the "interaction" of m supermatrices is presented and 
calculated explicitly in the cases of 2- and 3-matrix "interaction". The suggested 
technique allows us to calculate both the spectral correlations and the correlations 
of the eigenfunctions taken at different energies and in different space points. 
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I. INTRODUCTION 

A. Conventional and unconventional Random Matrix Theories 

The Random Matrix Theory (RMT) is a very useful mathematical formalism which allows 
to describe universal properties of complex quantum systems. Let us consider an ensemble of 
N x N Hermitian matrices, whose elements are independent Gaussian distributed random 
variables with a zero mean value and a position dependent variance: 

W = 0; <fl2> = i, (\Hi#\ 2 ) = ±&r{\i-j\)- (i) 



Here (. . .) denotes averaging over different realizations of RMs; the parameter (3 corre- 
sponds to the Wigner-Dyson symmetry classes: (3 = 1 for the Gaussian orthogonal ensem- 
ble (real matrices) and (3 = 2 for the Gaussian unitary ensemble (complex matrices). The 
function T and the parameter B determine various universality classes sharing the same 
global symmetries. 

A special case of the constant variance of the off-diagonal elements B 2 Ti\i — j\) = 1 
corresponds to the archetypal Wigner-Dyson RMT [lj. It has a great number of appli- 
cations starting from nuclear physics to quantum chaos to mesoscopic physics^. Ergodic 
wavefunctions and a level repulsion are essential features of the Wigner-Dyson RMT. 

Recently, unconventional RMTs characterized by decreasing function T have attracted 
a substantial interest (see, for instance, detailed introductions in references {4,5.6]). This 
interest is stimulated, in particular, by a possibility to explore the properties of localized and 
critical disordered systems. For example, if the off-diagonal matrix elements are essentially 
non-zero only inside a band centered at the main diagonal and decay exponentially fast 
to zero outside the band, all eigenf unctions are exponentially localized. This banded RMT 
describes the physics of a quasi-one-dimensional disordered wire^. If T decays only as a 
power-law outside the band^, T ~ lj\i — j\ 2a , the eigenfunctions are power-law localized for 
a > 1. The Wigner-Dyson universality class is approached for a — > 0. Thus, the power- 
law banded RMT (PLBRM) can interpolate between the Wigner-Dyson statistics and the 
Poisson statistics of localized system. The special case a = 1 corresponds to a critical 
behavior similar to that found at the point of the Anderson metal-insulator transitions*^. 
The function T C rit for the one-parametric family of the critical PLBRM can be defined as 
follows: J- critic — j\) — V (-S 2 + (i — j) 2 ) ■ The eigenfunctions of this critical model remain 
multifractal at any B ranging from the weak multifractality at the large band- width B ^> 1, 
to the strong multifractality for the almost diagonal RMT B <C 1. 



B. From large to small bandwidth RMT: cr-model versus virial expansion 

The considerable progress in the banded and the power-law banded RMT has become 
possible due to the mapping^ onto the nonlinear supersymmetric a-model 11 , which is a 
powerful filed-theoretical description of various averaged correlation functions. However, 
such a mapping is only justified in the large bandwidth limit. This limitation comes from 



the saddle-point approximation which is a crucial step in the derivation of the cr-model. 
Physically it corresponds to the diffusive approximation which implies in particular that 
the smooth envelope of a typical eigenfunction changes slowly on the scale of the mean free 
path. This approximation fails in the opposite limit of the small bandwidth, including the 
case of the almost diagonal RMTs where the bandwidth shrinks to zero and the off-diagonal 
matrix elements are assumed to be parametrically smaller than the diagonal ones, B <C 1. 

Let us represent the almost diagonal matrix as a sum of the diagonal part and a matrix 
of the small off-diagonal elements H = + V , V ~ B <C 1. The diagonal matrix 
represents "non-interacting" energy levels or localized eigenstates. The presence of the small 
off-diagonal matrix V leads to a weak "interaction" between different localized states. In 
order to calculate correlation functions, one can perform an expansion in the number of 
interacting localized states^. The small parameter B is the control parameter of this 
procedure. This method was called^^ "a virial expansion" (VE) by analogy with the 
expansion of thermodynamic functions of a dilute system in powers of density with the 
m-th virial coefficient being governed by collisions of m + 1 particles. 

These ideas were initially implemented in a semi-empirical real-space renormalization 
group approach which has been applied for critical systems with long-range interactions^ 
and for a quantum Kepler problem^. The real-space renormalization group was also used to 
study the critical almost diagonal PLBRM^: the spectral correlations and the scaling prop- 
erties of the eigenfunctions were investigated by considering a resonant interaction of two 
energy levels. The renormalization group approach has however two serious disadvantages: 
(i) it does not allow for a rigorous control of omitted contributions; (ii) due to technical dif- 
ficulties, it is almost impossible to go beyond the leading term, thus the resonant interaction 
of three and more levels remains inaccessible in this framework. 

A rigorous counterpart of the renormalization group approach has been suggested in re- 
cent works 0J^Q . It deals with a regular VE generated with the help of the Trotter formula 
(TVE)^. Based on the classification of perturbation series by a number of the interacting 
levels involved, TVE allows one to study the density of states and spectral correlations of the 
almost diagonal RMTs. The accuracy of the TVE is always controllable resulting in the rig- 
orous perturbation theory. The second disadvantage of the renormalization group has been 
also partly overcome: the TVE allows to go beyond the leading term by considering the in- 
teraction of two and three levels. In this way the first and the second virial coefficients have 



been calculated for a generic model of the almost diagonal RMTs. The general formulas 
have been applied 4 -^ to the different models of almost diagonal RMTs, including criti- 
cal PLBRMs^, the unitary Moshe-Neuberger-Shapiro model 1 -, and the Rosenzweig-Porter 
model 16 in the regime of crossover. 

The TVE involves a complicated combinatorial part of intermediate calculations. The 
combinatorial problem appeared there resembles a coloring problem of closed graph edges. If 
one considers the density of states^, the coloring is similar to the well studied problem of the 
graph theory^ whose solution is known 1 ^. However a study of spectral correlations requires 
to resolve much more complicated problem of simultaneous coloring of several graphs 4 . A 
complexity of the combinatorial calculations grows tremendously with increasing the num- 
ber of the interacting energy levels. Therefore the TVE can be used practically only for 
calculation of the first and the second virial coefficients. 

C. Virial expansion from the field theoretical representation 

In the present work we formulate a super •symmetric field-theoretical representation for the 
VE of different correlation functions of the generic model ([T]). The method of the super- 
symmetry allows us to perform an averaging over RM ensemble for an arbitrary function 
— The super symmetric VE is controlled by a small parameter B C 1. The 
virial coefficients are straightforwardly derived in a general form in terms of the integrals 
over supermatrices. It is important that no combinatorics appear in the intermediate field- 
theoretical calculations. In this framework, the interaction of m energy levels is described 
by an integral containing only m independent supermatrices associated with m different 
sites in the real space. In order to calculate the integrals over supermatrices explicitly we 
employ a parametrization introduced recently in Ref . [201] . 

The suggested super symmetric field theory (SuSyFT) can be equally applied both to the 
spectral correlations and to the correlations of the eigenfunctions taken at different energies 
and in different space points. We emphasize that this approach is the unique analytic tool 
to describe the wavefunctions correlations for the almost diagonal RMTs. Neither the real 
space renormalization group nor the TVE are capable to do this. We would like to mention 
also that SuSyFT might be useful for non-perturbative calculations as well, however this 
issue is beyond the scope of the present work. 



The paper is organized as follows: we present the SuSyFT and give main definitions in the 
Section [III The basic ideas and parameters of the VE are explained in the Section II II I with 
a reference to SuSyFT. A general integral expression for m-th virial coefficient governed by 
the interaction of m supermatrices is presented in the Section IIVI In this section we also 
discuss in detail the validity of a saddle-point integration over massive degrees of freedom 
and the applicability of the VE. We exemplify the integral calculations for the cases of 2- 
and 3-matrix interaction in the Section [V] and discuss advantages of the method and possible 
further applications in Conclusions. 



II. MAIN DEFINITIONS 



Let us introduce the retarded and advanced Green's functions 

G R/A (E) — 1 ; (2) 

E - H ± tO 

H is a Hermitian RM of large size N ^> 1 . It has independent matrix elements and belongs 
to the Gaussian ensemble described by Eq.([T]). Without loss of generality, we consider the 
case of GUE, (3 = 2. A generalization to other symmetry classes is straightforward. The 
Green's functions are an efficient tool to study different correlation functions. For example, 
the expression for the averaged density of states in terms of the Green's functions reads 

1 N 1 
(p(E)) = „ - en)) = ^ (Tr (6 A (E)) ) , (3) 

n=l 

where e n are eigenvalues of the random matrix H and (...) denotes averaging over the 
ensemble of random matrices. The inverse density of states taken at the band center E = 
governs the mean level spacing of RMT: A = 1/A(p(0)). The mean level spacing of almost 
diagonal unitary RMT is 5 

A 



~ (v^ + O^ 2 ))^- 1 

B<1 



The two-point correlation functions can be expressed by means of the quantity Q pq {uj): 

g m (u) = - / dE G* p {E + u/2)Gf q (E - u/2) . (4) 
For example, the expressions for the averaged two-level correlation function 

/oo 
dE((p(E + u/2)p(E-u/2))) (5) 
■oo 



is given by 

2ti 2 N 



2 N 



p,q=l 



where ((ab)) = (ab) — (a)(b); while the averaged correlator of two eigenfunctions taken at 
different energies at different space points p and q 

C 2 (iu,p,q) = A / dE((J2 S(E + u/2-e n )5(E~uj/2-e m )\A n (p)\ 2 \Ajq)\ 2 )) (7) 

m,n=l 

reads 

C 2 (u,p,q) = ^K((g pq (u J ))). (8) 
Thus, we have to calculate ((Q pq (uj) )) to explore the two-point correlation functions. The 



ensemble averaging can be performed with the help of SuSyFT lllJl9l |. To this end, we 
introduce iV supervectors 



/ ; a = l,2,...iV. (9) 




Here, ^ R ' A (a) = (s R ^ A (a),x R ^ A (ce)) T are supervectors in retarded-advanced sectors. They 
consist of commuting s R l A (a) and anti-commuting Grassmann variables x R ^ A ( a )- The 
direct product of \l/(a)-vectors with the conjugated vectors constitutes N Q-matrices of the 
size 4x4. The block structure of the Q-matrices in the retarded-advanced notation is: 

Q a = tf(a ®¥(a = KJ K \ K KJJ .],K=\ .(10) 

The ensemble averaged Q pq can be written as follows 



N \ / N 



(GvM)= A f A E j V{Q}K p A q (jQe^j [II e Sp[QM j ; (11) 

where T>{Q} = Yla=i ^{Qa} is the measure of integration over the supermatrices Q a (see 
Appendix |A] for the details)^. The factors 1Z P and A q break the supersymmetry between 
the commuting and anti-commuting variables in the retarded/advanced sectors: 

n p = ( X R (p)Y X R {P) , A q = ( X A (q)Yx A (q) ■ (12) 



In Eq. liTTj) we have separated out two parts of the action: S corresponds to the diagonal 
part of RM and depends on a single supermatrix 



So[Q a ] = Str^-aQl + i^E + ^Q^ , A=r°j , (13) 



RA 

a = \{Hl a ) = \, n^ui + iO. (14) 

The second part of the action S p is proportional to the variance of the off-diagonal elements 
and contains a product of two supermatrices 

Sp[Qk,Qm] = -hmSte{QkQrn} , k^lTl, (15) 

b km = l -(\H km \ 2 ) = - A B 2 F(\k-m\). (16) 

In the other worlds, S p describes an interaction of different Q-matrices. Note that 
Egs-f llllfTBTrTBl) are exact. We have used Q-matrices to compactify the notation but we 
could equally write the expression for (Q pq (oo) ) in terms of the integrals over ^f- vectors. 

The standard cr-model derivation includes (i) introducing auxiliary supermatrices 
that allow to decouple the quartic in (i.e., quadratic in Q a ) terms by the Hubbard- 
Stratonovich transformation; (ii) Gaussian integration over the supervectors After 
these two steps, one has to employ the saddle-point approximation for the integrals over 
. As we have already discussed in the Introduction, this procedure is justified only for 
the large band width RMT, B ^> 1. Here in contrast we will consider the case 

B<1, (17) 

where the standard a-model fails. Therefore instead of the Hubbard-Stratonovich trans- 
formation with further standard steps, we will use the method of the virial expansion and 
employ a parametrization of the Q-matrices suggested in the reference 0]. Details of the 
parametrization are given in Appendix We emphasize that integration manifolds of the 
<7-model and of the VE are quite different it is given by x u{i)xu(i) ^ or ^ e f° rmer 

and by jj^Q§n\ x K. for the latter—. In particular, a saddle-point approximation in the 
SuSyFT suggested will lead to the linear constraint Str(Q) =0 for a large scale theory 
(see Sect. HVT) whereas the saddle-point manifold of the standard diffusive a-model requires 
the additional nonlinear constraint (Q^) 2 = 1. Notice however that the nonlinearity of 
SuSyFT follows already from the definition of the matrix Q = \l/ <g> ^ . 




FIG. 1: Scheme of the SuSyFT virial expansion: (a) Supermatrices with different indices (as- 
sociated with different sites in the real space) do not interact with each other in the case of the 
diagonal RMT possessing the localized wavef unctions; (b) Pair interaction of of the supermatri- 
ces via (small) off-diagonal elements of RMs, which governs the first virial coefficient Q^ 2 \ see 
Eqs. (l40|4ip ; (c) The interaction of three supermatrices, which governs the second virial coefficient 
a (3) , see Eas. fjTOTD . 

III. BASIC CONCEPT OF THE VIRIAL EXPANSION 

The two-fold product of exponentials in Eq.f fTll can be expanded in power series 

(N \ oo 1 / N \ k oo 

JJ e s P [QM \ = J2 i J2 (-2& m ,nStr[g m Q n ]) ee V D + £ V<™> ; V D = 1 . 
n^m / k=0 \m>n=l / m=2 

(18) 



We have rearranged this series by separating out terms y( m ) which contain a given number 
m > 2 of different Q-matrices. It is easy to show that 

N 

y(2) _ ^ V^ a , V^a = e ~ 2ba ^ Str[Qa ^ Qa ^ - 1 ; (19) 

ai>(*2=l 
N 

V (3) = £ VgU, (20) 

ai>02>c«3=l 

y(3) _ y(2) y(2) v (2) , y(2) y(2) , y(2) y(2) , y(2) y(2) . 

(2) 

and so on for the higher terms. The function Vaia 2 m VE is a counterpart of the Mayer's 
function used in the theory of imperfect gases 2 ^. 

This expansion can be now substituted into the expression for ( Q pq {uo) ) : 

N 



JV+1 poo 



A 



J°°dE Jv{Q}TZ p A q ( n eSolQa] ) (V D + V (2) +V (3) + ...). (21) 



The first term V D = 1 corresponds to the diagonal RMT with noninteracting localized 
eigenstates. This is the starting point for VE. The wavefunctions of the diagonal RMs are 
completely localized at different sites having no overlap with the other sites and the energy 
levels are uncorrelated: 

0, (22) 



7T / 2 I ^ 2 

n 

A 2 



«0 = i I We-™ + %\ -/) (c-i-fi +/)):(23) 



CO 

2 



R^oo = ^—((g^(u))))=A{5(u)+0(l))=8(s) + 0(l/N) ] (24) 



^- , erfi(^) = -^= / e^dt . 





This is reflected by the structure of the first term: the Q-matrices are decoupled and all 
superintegrals factorize. We can say that the matrix Q a is associated with the site a and 
the supermatrices at different sites do not interact with each other if we put S p = (see 
FigOtti). 

Each term y( m ) is the sum of the exponentials containing m supermatrices. The super- 
matrices are linked by the small off-diagonal elements of RMs (see FigJTb,c). We refer to 
these links as "the interaction of the supermatrices" . Obviously, they reflect the interaction 
between the localized eigenstates of the diagonal part of RMT. 



Let us discuss general properties of the summands Vj^V, {a} = an, a 2 , ■ ■ ■ , a m . The set 
{a} must include external indices (p ^ q and p = q in the off-diagonal case and in the 
diagonal one, respectively); otherwise the contribution is canceled by subtracting the 

decoupled term, see the definition of ((. . .)). The N — m supermatrices, whose indices 
do not belong to the set {a}, are included neither in y( m ) nor in the symmetry breaking 
factor TZpA q and, thus, are containing only in "non-interacting" part of the action So- The 
integrals over these supermatrices are equal to unity due to the supersymmetry. 

Let us introduce several definitions. We will call supermatrices Qp 1 and Q@ 2 entering into 
the expression for vte; connected, if and only if there is a sequence of indices {71, . . . , , 7„} 
such that all elements &/3 l7l , & 7l72 , . . . , & 7n _ l7 „, b ln p 2 are contained in vfe . We refer to (i) a 
subset {a'} C {a} as to connected subset if all independent matrices Q{ a >} are connected to 
each other; (ii) two non-intersecting subsets {a^} C {a}, {cr 2 )} C {a}, {d^} f]{a^} = 
as to disconnected subsets if none of the matrix of the first subset Qj a m\ is connected to 
any matrix of the second one Q{ a {2)y If one or two external indices belong to a connected 
subset (p E {a'} or g G {«'} or p, q G {a'} ) then all matrices of this subset are connected 
to the supermatrices with external indices. These connections break the supersymmetry for 
all supermatrices Q{ a '}- 

If m = 2,3 then there are no disconnected subsets: all matrices entering Vj™j are 
connected to each other and hence all of them are connected to the supermatrices with 
external indices and the supersymmetry is broken for all of them. 

Starting from m > 4, the disconnected subsets may appear and there are containing 
supermatrices, which are not connected to any supermatrix with an external index. These 
are counterparts of the vacuum diagrams in the standard interaction representation^ and 
they do not contribute to the correlation functions. For example: (i) if {a*- 1,2 **} are two 
disconnected subsets and p,q E {c^ 1 -*} then the supersymmetry is unbroken for all matrices 
{Q| Q (2)}}, the integrals over the supermatrices with unbroken supersymmetry yield trivial 
boundary contributions governed by Q{ a (2)y = and the corresponding Mayer's functions 
are zero; (ii) if {c^ 1,2 )} are two disconnected subsets, p 7^ q and p E {c^ 1 -*}, q E {a^} then 
the contribution V/™? is canceled by subtracting the decoupled term (see the definition of 

((•••)))• 

To summarize the above discussion, each term V*™- 1 effectively consists only of irreducible 
parts where there are no disconnected matrices. The supersymmetry is broken for all m 



matrices (i) by the factor 1Z A for the supermatrices Q p , q and (ii) by the inter-matrix links 
for the rest supermatrices. To formulate the VE, we focus just on these supermatrices 
Qa x , Qa 2 ■ ■ ■ Qa m with broken supersymmetry. Let us scale them by the constant B [24j : 



Q a . = BQ ajl l<j<m. 



(25) 



After this scaling, the small parameter B is eliminated from y( m ) but appears in the single- 
matrix part So in three different ratios as l/B 2 , E/B, and Vt/B: 



So 



Qc 



B 



Str 



(26) 



The first ratio is large, 1/B 2 ^> 1, and, therefore, the exponentials 



exp < S 



B 



oc exp 



4# 



Str 



Q: 



< B 2 . We can draw 



suppress the volume of integration over Q a . by constraint Str Q c 
a conclusion that the larger number m of independent Q-matrices the smaller is the con- 
tribution Q^ 1 of m-matrix term V^ m ^ to the correlation function. This statement would 
be unquestionable if one integrates only over commuting variables. In the case of superinte- 
grals it is more subtle and requires an additional discussion. An obvious counterexample to 
our estimate is the case of a partition function, for which the supersymmetry is unbroken, 
1Z = A = 1, and all superintegrals yield unity regardless of the apparent small phase volume 
of the integration. This happens due to the anomalous contributions^ to the superintegrals. 
The irreducible terms vfe can also contain the anomalous parts. The detailed analysis of 



the correlation function (see Section IV CI) shows however that, due to the broken supersym- 
metry, the latter anomalies do not change the powers of the small parameter B and the 
following estimate holds true: 

Mm+l) 

ypg ocB, m>2. (27) 



r (m) 

This justifies the supersymmetric virial expansion for the almost diagonal RMT: the ex- 
pansion of the correlation function in the number of the interacting supermatrices, Eq. 
yields a regular perturbation theory in powers of B. 



IV. GENERAL EXPRESSION FOR THE VIRIAL COEFFICIENTS 



A. Saddle-point integration 

Before presenting an expression for the arbitrary number of the interacting matrices, let 
us analyze the 2- matrix term {{Qpq )) in more detail. As we have already mentioned the 
supermatrices contained in V*- 2 -* must be coupled to the factors 1Z P and A PtQ . Therefore, 
we obtain the following expression for the off-diagonal p 7^ q and the diagonal p = q parts 
of «<$>»: 

((C)) = -X ( 6 ^tr[Q p + Q q }) K p A q V$) , (28) 
? N 

((O) = "X E (HMQp + Q*])Kp>ApV<$) • (29) 

n=l,n^=p 

We have introduced the averaging over the supermatrices: 

I V{Q a }{. . OtAPM ; £ [Q Q ] = S [Q a ] =-^^ + i %r[AQ a ]. (30) 
J e=o 4 2 

The ^-function in Eqs.( l28ll2"9"]) resulted from the integral over E, see Eq. (ll3|l21l) . 

(2) 

Formulas for Q pq exactly describe the contribution of the two-level interaction to the cor- 
relation function Q pq . However, contains a product of two supermatrices, cf. Eq. (fT9l) . 
which entangles the integration variables in a nontrivial way (see Eq. ( IAIOI) ). We will calcu- 
late the integral over the variables R a (see Eq. (IA8l) ) in the saddle-point approximation. To 
explain this step, it is convenient to scale the supermatrices by ^Jh p ~ OL 

Qp \J^pa Qp j Qa \J^pa Qa j 

consider the integration over the scaled variables R a and S a in Eqs. fl28ll29p and perform 
the approximate integration over R: 

R 2 B 2 

00 JJ\Rp, a \ \ \fOpa J s p - n p D a — rt Q 

^ (31) 

^^(jK W ^^^L ei *^ +Sa)v S(°' 5 ) +°(v / Q + o (|) 

£ a = Str[Q a ], 5„ = Str[AQ a ]. 



Here a = n,p for the diagonal- and off-diagonal parts of Q^ 2 \ respectively. We have 
accounted for the 5-function in Eqs. ( I28f29l) and denoted the typical values of S and R, 
at which the integrals converge, by S t and R t accordingly The value of R t is fixed by 
the Gaussian exponentials in Eq. f[3"Tj) : R t ~ \fbp~ a ■ The integrals over S converge due to 



<> 



-s. 



the exponentials e 2 ; cf. detailed calculations in the Section |V] Therefore, we can 



estimate the second characteristic scale as S t ~ yhpa/oj |25|] and obtain Rt/S t ~ uj. The 
saddle-point integration over R makes sense only if the corrections in the right-hand part 
of Eq.(EH) are small^: 

max | a/?V, ~ max j \[hp~ a , w| <C 1 . (32) 

Thus, we have to restrict ourselves to the region uj <C 1 where the density of states of the 
almost diagonal RMTs is close to the constant^. The ultrahigh frequencies uj > 1 cannot 
be considered within the saddle-point integration over R and they are beyond the scope of 
the present paper. We can return to the unsealed matrices Q and arrive at the following 
equation for : 

«0 = -2^(*(Str[Q p ]) *(Str[Q ff ]) n p A q V$) +Sg^ g , (33) 

/o \3/2 N 

= - [ -^— 53(*(Str[Q p ]) 5(Str[g n ]) 1l p A p V$) +6g$. (34) 

n^p 

To calculate the leading terms, we have effectively replaced the exponentials exp (— R 2 /4) 
in the integrand of (I28H30P by the 5-functions of R: 

e^f-istr[Q; + Qy)*^ (35) 
V2^5(R a ) 6(Rp) = V2tt £(Str[Q a ]) 5(Str[Q p ]) . 



Eqs. ( }31f32j) implies in the case of the off-diagonal correlation function: 

-^~max{Sv/^(b-g|),^} . (36) 
A similar estimate for the diagonal correlation function is more subtle since the derivation 

(2) 

of Opp involves the summation over the auxiliary index n. Let us assume that this sum in 
the leading part of Qpp as well as in the correction SQpp converges at a characteristic scale 
|n — p\ ~ X c , then we can expect that 



r2) max{By/F(X c ), u} . (37) 



The value of X c is, of course, model dependent and varies for different RMTs. 

The generalization of Eqs. (l33|l34l) for an arbitrary number m of the interacting Q- 
matrices reads 

«0 * 2 ~K h2V T' 1 £ (s(MQ P ])HMQ q })K P A x (38) 
II^(Str[g Q ,]) ) VM ia2 ... Qm _ \ l<j<m-2; 



im-l ^ 



«<$°» - f ( £ ( 5 ( Str ^D >< ( 39 ) 

II^(Str[g a ,])J VS ia2 ... Qm _\ l<j<m-l. 

' ' QpiQq Qa.2 •■■Qa m _] 

The summation is performed over ordered indices 1 < oi\ < ct<i < a% < . . . < N excluding 
the external fixed indices p and q. It is easy to show (see sect ion IV CI) that these expressions 
have the following functional dependence on parameters Q and B: 

((g {m) )) = B m - 2 g {m) (n/B) , (40) 

with function depending only on the ratio Vt/B. Thus one can write down VE of the 
correlation function as a functional series in powers of B: 

((g pq (B, n/B))} ~ Q° + ® m ~ 2 &£\n/B) . (41) 

m>2 

The functions Qpq are the virial coefficients. Each coefficient Q( m \Q/B) is governed by 
the interaction of m supermatrices corresponding to the interaction of m localized states. 
The first term Q D is related to uncorrelated statistics of the diagonal part of the almost 
diagonal RMTs. 

By analogy with the estimates (13611371) . Eqs. (l38|l39l) describe an arbitrary virial coefficient 
g( m ) with the following accuracy: 

-^~max{i3^(b-g|),^}, -^-maxjtfv/^),^}. (42) 
It is convenient to represent the corrections schematically as a sum of two terms: 



where S u gff /gff ~ u, hQ^JQ^ ~ By/F(\p-q\), and 8 B $? ~ B^VQ- 



B. Validity of the virial expansion and the large scale limit 



All details of the evaluation of C^ 2 ' 3 - 1 and are presented in the next section. Here, we 
would like to discuss the validity and the applicability of VE ( 138MB . 

Firstly we note that VE (1411) is a functional series. Its successive terms decrease with 
increasing the number of the interacting supermatrices only if the absolute value of the virial 
coefficients is bounded for the arbitrary ratio Q/B. This condition determines a convergence 
of VE but it cannot be checked until the RMT model is specified. In particular, this condition 
is violated for RMTs with almost ergodic wavefunctions. 

Secondly we should recall that, calculating (Q/B) by the saddle-point approximation 
in R- variables, we have neglected corrections S^Qp^ which are of the order of 0(u). This 
means that, for a given u the summation over m in VE described by Eqs.( l38H4T]) must be 
stopped at 

l + log(u;)/logOB). 

For instance, if B < lo <C 1 then 5 U Q^ > and one may take into account only the 
interaction of 2 supermatrices neglecting all higher terms. The next term of VE governed 
by the interaction of 3 supermatrices may be taken into consideration only for the smaller 
energy cu <^ B when d^G^ <C G^\ etc. On the other hand, the neglected dependence of 
the virial coefficients on u often results from the energy dependence of the density of state 
and does not influence universal properties of the correlation functions. 

Finally let us compare 5 B g^ with the successive term of VE £( m+1 ). For the sake of 
simplicity we compare SsGplq with Gpl q though the same analysis applies to the diagonal 
virial coefficients with p = q and for the arbitrary m. Without loss of generality we put 
= 1. If \p — q\ ~ 1 then $BGpl q /G^l q ~ B and there is no way to get the scale 
separation: 5qQ^ and Q^ 1 are of the same order and, again, one may consider only the 
two matrix interaction regardless of the energy smallness. 

Now, we will show that 5bG^ is parametrically smaller than the higher terms of VE in 
the large scale limit. The large scale limit means that we consider only those correlation 
functions which are not sensitive to the details of T at small distances and governed by the 
behavior of this function at large distances X c , at which T is sufficiently small T{X^ <C 1. 
It means, in particular, that we assume: (i) \p — q\ > X c in the case of G^lq] (ii) the main 
contribution to the sum over the auxiliary index a± originates from \cti — p\ > X c and 



\oti — q\ ~ X c in the case of Gp2 q - We remind that the sum over ct\ is due to presence of 
the third supermatrix, see Eg. (I38j) . Although the second assumption is not applicable for 
arbitrary function T , it allows to study a wide class of almost diagonal RMTs. We arrive 
at the following estimate in the framework of the large scale limit: 

;(2) 

BT{\p-q\)<BT{X c )<^B. 



Qp+q 



The ratio Qpl q /Qp^ q requires a separate consideration: the presence of the third super- 
matrix in the expression for G^ q results in the additional factors B^f T{\p — aij) or 
ByJ~T(\q — a\) , and simultaneously requires an additional summation over the index ol\, 
see detailed calculations in Sect lVBl below. This summation is crucial: we have assumed 
that it converges at \p — m\ ~ \q — m\ ~ X. , but the large phase volume of summation can 



compensate the smallness of a/^(X c ), cf. [6]. If this is the case we obtain 

(2) 

Thus, if the assumptions of the large scale limit hold true then 5sQ p ^ q is parametrically 
smaller than the next term of the virial expansion. In particular if the characteristic scale 
X c depends on iV and the function T decreases in such a way that limjy^oo T{X^) = 
then the ratio 8sQ^l q l 'Qp\f q asymptotically goes to zero. This scale separation justifies VE 
in the large scale limit in many cases. One can check, for example, that it is correct for 
the spectral statistics of the critical almost diagonal PLBRMs, where X c cx BN and the 



relevant energy range reads uo < BA {4]. 

We would like to mention that the large scale limit considered here is analogous to the 
diffusive approximation of the standard cr-model. In the latter approximation the spatial 
scales large compared to the mean free path are assumed to be the only relevant ones. In the 
same time the saddle-point approximation justified in the large scale limit in our approach 
results in the linear constraint on the Q-matrix 

Str[Q] = 0, (43) 



while the saddle-point manifold of the standard diffusive cr-model is defined by the additional 
nonlinear constraint (Q^) 2 = 1. 



C. Correlation function at the band center 



The expressions for the correlation functions given by Eq. (138|) and Eq. ()39p were derived 
after integration over E. However in some applications it is more convenient to consider 
a correlation function at a fixed energy. The aim of this section is to discuss briefly the 
correlation functions at the band center E = 0. 

We define the two-point correlation function at the band center as 

g pq (u, E = 0) = ££( W /2)G£(^/2) . (44) 
Repeating all the steps leading to the results Eq. (I38|l39j) . we obtain: 

({G™(u,E = 0)))^^({g£\u;))). (45) 

The ratio ^frajN can be referred to as "the unfolding factor"—. 

Let us also note that an average of the product of two retarded (advanced) Green's 
functions can be neglected again for the following reason. If we look at the parametrization 
of Q in this case (Appendix [B]) , then we notice that variables R and S change their roles. 
For this reason the large scale approximation Str[Q] « implies now \r, \r> ~ 0. Thus the 
volume of the integration in this case becomes parametrically small at B <C 1 and w< 1. 

Below, we will analyze only the correlation functions averaged over E. 

V. THE CASES OF 2- AND 3-MATRIX APPROXIMATION 

In this section we present detailed calculation of the contributions and governed 
by the interaction of 2 and 3 supermatrices respectively. It is more convenient to expand 
V (2 ' 3) in the series of powers of Q-matrices and then integrate over Q term by term. This step 
is not essential for the approximation of two interacting matrices which can be worked out 
directly from Eqs.( l33|34j) . However, it is more convenient for the 3-matrix approximation, 
since it allows one to unify the calculations for the different number of the interacting 
supermatrices and to avoid an explicit derivation of anomalous terms in the superintegrals. 
The disadvantage of this route is that the series obtained converge only asymptotically and 
one has to Fourier transform them in order to analyze the answer in the time-domain^. 



A. The case of 2 interacting supermatrices 



We start with calculating using approximate formulas (13311341) . The power series for 

l/ 2) reads: ^ 

y(2) _ (~2&pnStr[Q p Q n ]) 

fc=i ^' 

We use the phase and two non-compact variables R and S to parameterize the boson- 
boson sector of each supermatrix (see the corresponding definitions in Appendix EJ. The 
integration measure in Eqs. (133f34p takes the form: 

J P{Q}*(Str[Q])(. . .) ^2j°°dR5(R) J™ ^ i- jf% | d{r^> A }(. . .) . (47) 



and the expressions for 7^pAj,p read: 

1 



1Z A 



4 

Integrating over i?-variables we obtain 

(2tt) 3 / 2 ^ rr°° dS, 



A S P S q (v*rVr) p (v*aVa) q ] (48) 
l S 2 p {vrVrV*aVa) p . (49) 



((C)) - E jr^/ q{Sp+Sq) J WnVRVlVAUx 

x to) p (^) g J£ ^- 2 \^ f , (50) 



..( * * \ /"/" d0 p n (—2b pn S pn ) rri\ 
X(V R VRV A VA} P jj o j^y 2 ^ • (51) 



The expression for <S M = StrKLQJ oc S'pS'o is given in Appendix [Dj Eqs.( 1D3IID5l) . 
The integrals of S£ over the phases are calculated in the same Appendix, Eqs. (ID2IID4l) . 



pq 

The integrals over the S'-variables are regularized at the upper limit by the imaginary part 

(2) 

of f2 and converge at the lower limit for all k in the case of G p -i q and for k > 2 in the 
case of Qpp . The term with k — 1 in the diagonal part Qpp is special: it is governed 
by an anomaly, i.e., an uncertainty x oo with zero resulting from the integrals over the 
Grassmann variables (vrVrVaVa) n and infinity due to the divergence at the lower limit of 
integration over the commuting variable S n . This uncertainty can be resolved either in a 
standard way-ii or, equally, one can calculate the integrals in the diagonal part for k > 2 



and then perform an analytic continuation for k = 1. The result of the integration over all 
variables can be written as follows: 



«0 

((Gg)) 



(2tt 



,3/2 



A 

(2n) 3 / 2 



E 



A 



k=l 

N oo 

EE 

n^p k=l 



(2b m \ k Y{2k-l) 



T(k) 



(k-1) 



(2b pn \ k Y(2k-l) 
V{k) 



\9? 



k. 



(52) 
(53) 



The correlation functions R<± and can be calculated from the real part of Q pq (see 
Eqs.( l6|8il ). Taking the real part by substituting u instead of Q in Eqs. ( l52|53l) one obtains 
the asymptotic series in the energy representation. However it is more convenient to consider 
the time representation by performing the Fourier transform of the real part 



GpS) = ^J dcje-^(g pg (cj) + c.c.), 



obtaining 



k=l 

N oo 



(^g } (*)» 



A 2 |t| 
7r(2vr) 3 / 2 vv 

2-^2-^ (k-l)\ 2k -1' 

n^p k=l v > 



, 26 pn t 2 ) k 



A 2 lt| 



Now the summation over k can be done explicitly 

-7T 5 / 2 



r 5/2 



\/26pg|t 



7T 



erf ( V26 M |t| 



5/2 ^ r /— 

«^(*)» * -V2^J]V2V v^l^- 2 ^ + ^erf (v^|t 



nj^p 



(54) 
(55) 

(56) 
(57) 



Here erf (2) = J* e~ t2 dt. Our theory can be verified by comparison with the results 
of TVE. To this end we calculate the form factor, which is the Fourier transform of the 
two- level correlation function R2: 

A 2 



N 



K(t) 



2n 2 N 



(58) 



p,q=l 



and insert in this formula Eqs. (I56|I57|) . This gives the form factor in the approximation of 
two interacting levels: 



N 



N\t\ 



E X (\p-Q\) e ~ 



x(\p-q\) 



p,q=l 



2n 



N 



-2—— > x(m)e 



(m) 



(59) 



m=l 



x(b-ffl) = 26 w t 2 = -(i3t) 2 ^(|j9-g|); 



which coincides with the expression for K^ 2 '(t) obtained by TVE 4 . 

This comparison of TVE and the theory based on SuSyFT clearly demonstrates that 
SuSyFT is capable to give much more detailed information on the correlation functions. 
Namely, TVE deals with the form factor which is an integral quantity obtained after the 
summation of diagonal and off-diagonal parts of the correlation function Q over all spatial 
coordinates, while the correlation function Q at given spatial points can be derived only 
from SuSyFT. 

We can now return from Eqs. ( 15611571) written in the time-domain to the energy represen- 
tation of ft [gw] : 

»{(O w ))) : 

The power series Eqs. ( I52|[3HT) are asymptotic expansion of these formulas^,. Note that the 
summands in the right hand side of Eq. floTl) are peaked around the value 

(62) 




V^I 2B^jn\p-n\) 

see FigfSJ Thus, we can find the characteristic spatial scale X c , which yields the main 
contribution to the sum over n and determines Gpp(s), from the following estimate 

HX C ) ~ (|) 2 • (63) 

Estimate ensures the validity of the large scale limit for uj <C B <C 1 at the level of 
2 matrix approximation. Indeed if uj <C B then ^{X^ <C 1 , hence X c ^> 1, i.e., the 
diagonal correlator Qpp is governed by the large distances, and the correction Ss to 
the saddle-point integration is expected to be smaller than the higher terms of the VE. On 
the contrary, in the range B < uj 1 the characteristic scale is small, X c —>■ 1 and the 
higher terms of the VE can be of the same order as the omitted correction SbQ^ 2 \ see the 
Section IIVBI This means that we cannot use the saddle-point integration to go beyond the 
two-matrix approximation in the case B<w<l. 




FIG. 2: Summand in the right hand side of Eq. (|61|) : f(y) = — 1 — ^-e y2 [2y + -) erfi(y) 



Here the argument y denotes the parameter u/y/8b 



pn- 



B. The case of 3 interacting supermatrices 

The calculations of based on Eqs. ()38|l39j) are very similar to those described in the 
preceding section for Q^ 2 > . The power series for read 

oo oo oo oo 

vi± = { V + V I + V I + V I 

I ^ ^ lfe 3 =0 ^ \k 2 =0 ^ lfcx= 

fel,2,3=l ^1,2 = 1 «1,3=1 fo,3=l 

(— 2b pn Stx[Q p Q n \) 1 (— 2bp m StT[QpQ m ]) 2 (— 26 mn Str[Q m Q n ]) 



(64) 



hi 

We insert this series in Eqs. 

, 2 N oo 



k 2 \ k 3 l 
and integrate over /^-variables obtaining 



«0 



(An) 



V3A 
x (VrVr) p {vaVa) 
(4tt) 



E E 

{m^p,q} fci,2,3=0 
2tt 



SpSgS, 



e 2 



d{V*RVRV*AVA}p,q,m x (65) 



^p,q,m ( 2hpqSpq) ( 26p m iSp m ) ( 2bq m Sq m ) 



({Ggh) - 



« llt (2vr)3 

2 TV oo 



E E 



hi 



d{^?7ii^^}p,m,n x (66) 



X 



(VrVrVaVa). 



2n 



■'p,m,n 



2bpqSpq) ( 26p m iSp m ) ( 2,bq m S, 



qrn ) 



(2tt) 3 fci! 



m > n . 



We use the same trick with the analytical continuation from kj > 2 to kj = 1 to handle 
the anomalous terms. Note that we have put zero in the lower limit for the sum over ^1,2,3 
and combined all 4 contribution in Eq. (164j) together. This is possible since all terms with 
either ki — k 2 = or k x = k 3 = or k2 = k 3 = are equal to zero after the analytical 
continuation (see the results for below). The integrals over the phases are calculated 
in Appendix [Dj see Eqs. (lD6HD12p . These rather cumbersome expressions are substantially 
simplified after the integration over the Grassmann variables (see Appendix |E|) : 

J d {VRVRV*AVA}p,q,m(VRVR) p (vlVA) q JjJ ^Pl^ OW* (Sgmf 3 = 



J d{VRVRVAVA}p,q,m(v*RVRVAVA) p JJJ^ ^ Spm ^ = 

^ gk-L+k-z gki+k 3 gk 2 +k 3 

32 E{ku h > h) r(jfe! + h - 1) v{k x r(k 2 + k 3 - i) ; (67) 



where 



*(k k k) - r(fci-i/2)r(fc 2 -i/2)r(fe 3 -i/2) 

x {2k x k 2 k 3 - k x k 2 - k x k 3 - k 2 k 3 ) T(k l + k 2 + k 3 - 1) . 
After the integration over S'-variables, the power series for take the following form: 

<«&» = -i^S t E (^)" (^f >< W 

{m^p,g} fci, 2 ,3=0 V / V / \ / 

x E(k x , k 2 , k 3 ) (ki + k 2 -l) (ki + k 3 - 1) , 

<«£»> = - S§ t E (^)" >< (*» 

V {m,n^p} fci, 2 ,3=0 V / \ / \ / 

x S(fci, k 2 , k 3 ) (h + k 2 )(ki + k 2 - 1) . 

Following the procedure described in the preceding section, we Fourier transform the real 
part of Eqs. fl69|70l) and obtain: 

3 i TV oo 

«&))^^ E £ HM 2 )* 1 Hwf Hw 2 )* 3 x (7i) 



4v/3 (iA) 2 

{m^p,<j} fci,2,3=0 

T(2[k 1 + k 2 + k 3 ]-l) y 1 ni ' 



3 N oo 

{{g ™ m ^iMuh £ £ Hw 2 )* 1 Hw 2 )* 2 x (72) 

{m,n^p} fci,2,3=0 

r(2|fci + fc 2 + fc 3 J - i) 

The triple sums on the r.h.s. of Eqs. (j71l72j) can not be reduced to a product of simple sums. 
Therefore, the summation over ^1,2,3 is not trivial. To verify SuSyFT, we calculate the 
contribution of three interacting matrices to the form factor (158j) . The answer looks more 
compact if we at first symmetrize the expression for G^(t) with respect to fci.,2,3 and then 
turn to the ordered sum over 3 remaining indices: 

N oo 

KW{t)*-±-± J2 E HM 2 )" 2 >< (73) 



{m>n>p} fci,2,3=0 

4 fci+fc 2 +fc 3 -i n E(kx,k 2 ,k 3 



(hi + k 2 + k 3 - l)(ki + k 2 + k 3 - 3/2) 



r(2[k 1 + k 2 + k 3 ] - 1) 

Eg. (1731) coincides with the expression for K^\t) obtained by TVEi. 

We remind that all results of this section should be taken into account in the VE if 
uj <C B <C 1 and the assumptions of the large scale limit hold true, i.e., the sums over m 
and n converge at the large spatial scales. This is the case, for instance, for the spectral 
statistics of the critical PLBRMs where the relevant energy range is small uj < B A while 
the characteristic spatial scale, which governs the two level correlations in the framework of 
2- and 3-matrix approximation, is large X c ~ B/uj 3> 1 {(J. 



C. Verification of the supersymmetric VE 

Let us verify that the power series obtained for , m > 2, really obey the estimate 
(1271) . If we consider a term of the power series with given powers kj , integrate it over all 
i?-variables in the saddle-point approximation, scale S"-variables by Q, and perform the 
summation over all internal indices m ,n ... ^ p,q, then a simple power counting shows 
that the answer will be proportional to 

(cf. Eqs. f l52|T5HT) and fl69|70l) ). Here Q m in the numerator and Q 2 in the denominator result 
from the integration measure and from the factor 1ZA, respectively. We can rewrite this 



ratio as follows 

g\ 2(fci+fc 2 +...)-(m-2) 

Obviously, having performed the summation over kj, the answer for can be written 

as a product 

g {m \B, n) = B m ~ 2 g {m) , (74) 

which agrees with the estimate (12 7p and with the formula (|40p . 

We remind that the m-th virial coefficient depends on the parameter Vt/B, which 

can take on an arbitrary value and is not assumed to be either small or large. The successive 
terms of the VE decrease with increasing m only if the absolute value of the virial coefficients 
is bounded for the arbitrary ratio Q/B. 




VI. CONCLUSIONS 



In the present work we develop a supersymmetric field theoretical description of a Gaus- 
sian ensemble of the almost diagonal Hermitian Random Matrices. In this ensemble the 
off-diagonal matrix elements are assumed to be parametrically smaller than the diagonal 
ones: Ha ~ 1 , Hij/Ha ~8<1. We use the method of the supersymmetry to perform an 
ensemble averaging. The standard route of the derivation of the supersymmetric nonlinear 
cr-model can not be taken in this case, since the diffusion approximation fails. 

As an alternative to the supersymmetric a-model, we derive a virial expansion (VE) 
in the number of "interacting" supermatrices, which is controlled by the small parameter 
B . Each supermatrix can be related to a localized eigenstate of the diagonal part of RMs. 
Thus, the supermatrix interaction describes the interaction of the localized wavefunctions 
via the (small) off-diagonal elements of RMs. The principle idea of VE is similar to one used 
in VE based on the Trotter formula^. Nevertheless, the supersymmetric VE is much more 
powerful since it allows to study not only the spectral correlations but also the correlation 
of wavefunctions taken at different energies and in different space points. 

The application of supersymmetric VE becomes especially efficient in a situation, when 
(i) the relevant energy range is much smaller than the typical value of the diagonal elements 
of RMs, u < 1; and (ii) the large scale approximation can be used, see Sect II VI In this 
case the massive degrees of freedom are integrated out by the saddle-point approximation. 



This step is a counterpart of the saddle-point approximation used in the derivation of the 
nonlinear a-model. However, the saddle-point approximation in the VE requires only the 
linear constraint: Str[Q] = 0. 

One of the main result of the paper is the integral expression for the m-th term of VE, 
Eqs. (l38ll39l) . which is governed by the interaction of m supermatrices. The superintegrals in 
this formula completely circumvent a complicated combinatorial calculations in the theory 
based on the Trotter formula. We note in passing that in this way we manage to reduce 
the complicated problem of simultaneous coloring of edges of several graphs (which is along 
standing problem in the statistical physics and the applied mathematics) to the calculation 
of superintegrals. The superintegrals are calculated explicitly for the cases of 2- and 3- 
matrix interaction with the help of the parametrization suggested in Ref . [20I] . The results 
containing in Eqs. (15211531) . (16011611) and (I69|70p have been obtained for the first time. They 
have been derived for a generic ensemble of the almost diagonal RMs described by Eq. ([T]) in 
the case of the unitary symmetry class. We note that our approach can be easily generalized 
to the other symmetry classes. 

The virial expansion generates a regular perturbation theory in powers of B for a variety 
of the correlation functions in the different models of the almost diagonal RMs. Applications 
to certain RMT models will be presented elsewhere. The critical ensemble of RMs^ and 
the Moshe-Neuberger-Shapiro model 15 with the orthogonal symmetry are two examples of 
promising applications. The next important step would be derivation of non-perturbative 
results directly from the supersymmetric action in the large scale approximation. The non- 
perturbative solutions could particularly shed light on the following problem: under what 
circumstances an interaction between the localized states can lead to the criticality or to 
the derealization. This question is of fundamental importance in the theory of disordered^ 
and strongly-correlated disordered systems^. 
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APPENDIX A: PARAMETRIZATION OF MATRIX Q IN THE 
RETARDED- ADVANCED SECTOR 



Let us find a parametrization of the matrix Q defined as direct product of the supervector 
by the conjugated supervector 

^ t t ® *h ® K^ a \ 

Q = ^®^=\ f f , (Al) 

*=(::) ■— fc) ■*■(-)■ 

where indices i? and A are referring to the retarded and advanced sectors correspondingly. 

Matrix Qrr = ^r<S)^ r is the orthogonal projector on vector ^r and can be diagonalized 
by the unitary matrix Ur\ 



(A3) 

II — ln%,ria —n* n \ , II — hn%,rin n%, \ I X 2 -, C 

U R = 



2 II 

Vr 




, Drr = | 













where 77/j = Xr/sr, \ 2 r = \\^r\\ 2 = \sr\ 2 + X*rXr- In the similar way block Q AA can be 
diagonalized by the pseudounitary matrix U A (U a KU a = K): 



(A4) 

/ 1 + hn*AVA 77* \ , ll + ln* A n A —n\ \ /-A 2 , 

U A = 






f-A 2 




, D AA = | 








I 





»M \-hf A r)A) \ -Va l-h* A r) A 



with 77^4 = Xa/sa, A^ = — H^aII 2 = Is^l 2 — X*aXa- Moreover Ur,Ua diagonalize off-diagonal 
blocks of Q: 

Qar = UaDarU r 1 Qra = UrDraU a 1 
D AR = ] D RA =[ , (A5) 

00/ loo/ 



u 



(A6) 



where e*^ = s* r sa/{\ s r s a {)■ Thus the matrix Q can be parameterized as follows: 

Q = UDU- 1 

( X 2 R -e-*\ R \ A 0\ 

e**\ R \ A -\\ 
0/ 

The measure in this parametrization can be easily found by calculating the Jacobian 
(Berezenian) of the transformation (1A6I) and is equal to 

(A7) 



(Ur 


0^ 


I- -1 


( D RR D RA \ 








1° 


U A J 




[D AR D AA ) 



\ 



V {Q} = — \ R - A d$dr\ R dr\ R drf A dr\ A . 



71 \ R X A 

Using Eq. flA6[) and taking into account the rotational symmetry of the supertrace, we obtain: 



Str(Q) = R, Str(Q 2 ) = R 2 ; 
Str(AQ) = S ; S = \\ + \\ ; 



R = \ 2 R - X\ ; 



(A8) 
(A9) 

(A10) 



Str[QQ] = A^(l - a* R a R ) + \ A \ 2 A (1 + a* A a A ) 

-2cos9X R X R \ A X A ^1 - ^a R a R ^j ^1 + ^«>a 

where a R/A = r) R / A - fj R / A ; 6» = </>-</> + A;A = ^{f} R r] R - rf R fj R + - t?a?U) and tilde 
marks the variables of the matrix Q. The last identity ( 1A10I) follows from the well known 
property^ 2 - of the matrices U : 

U R \m)U R { m ) = U R {T} R -fi R )e^^-^\ 



U A (fj A )U A (r] A ) = U A (r} A -f) A )e ? 



(All) 
(A12) 



APPENDIX B: PARAMETRIZATION OF MATRIX Q IN THE 
RETARDED-RETARDED SECTOR 

In the retarded-retarded sector the matrix Q is defined similar to Eq. flAip but without 
matrix K: 

• y R <g> * r ® *L \ 

Q = ® M/t = ( R R , (Bl) 

*= | | : Vr,r> = | | (B2) 




where indices R and R! are referring to the retarded sectors of two different Green's functions. 
As a result one can diagonalize Q by transformation similar to Eq. flA6p : 

Q = UDU- 1 
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(B3) 



where Ur is defined in Eq. (IA3p and Ur> is obtained from Ur by replacing subscript R by R' 
everywhere. All parameters appearing in (IB3P are defined in the same way as before: 



VR 
Vr' 



Xr 

Sr' 

XR' 
sr> 



X 2 
A R 



\sr\ 2 + XrXr, ^ 



S R S R> 
SRSR' 



\2 I |2 , * 

X R > = \sr>\ + Xr'Xr' 



(B4) 



The measure in this parametrization remains the same as in Eq. (IA7j) . while the roles of 
variables R and S defined in Eq. (IA8p and Eq. (lA9p are interchanged now: 

(B5) 
(B6) 



Str(Q) = S , Str(Q 2 ) = S 2 ; S = X 2 R + X 2 R , , 



Str(AQ) = R ; R = X 



R 



X 2 

A R> ■ 



Finally the expression for Str[QQ] is again similar to Eq. flAlOp : 

Str[QQ] = X R X 2 R (1 - a* R a R ) + X RI X R ,(1 - a* R ,a R >) 

-2 cos 6 XrXrXr'Xr' (l - ^a* R a R j (l- ^a* R ,a RI 



(B7) 



where a R / R > = T]r/r> - fj R /R> ; 6 



A ; A = Wr^r - v*rVr - Vr>vr> + vhVR')- 



APPENDIX C: INTEGRALS FOR INTEGER POWERS OF sin(</>/2) 
One needs the following formulas to average a product of the supertraces over the phases: 

Hh , h , h) s _i_ j£ sin *. (£) (|) sm - (^) ^ = (02) 
T(h + 1/2) T(k 2 + 1/2) T(fc 3 + 1/2) T(h + k 2 + k 3 + l) 

vr 1 / 2 vrV 2 ttV 2 r(fc! + k 2 + i)r(A;i + fc 3 + i)ro 2 + fc 3 + 1) • 



We consider only integer powers of sines k and k\^ 3 . Eq. flCip can be found in standard 
mathematical tables^ while Eq. flC2j) can be proven by the induction over one of the expo- 
nents, for example, over k 3 . Clearly, jF(fc,0,0) = J-{k) ; ^(ki, k 2 , 0) = T{k\)T(k-i). This 
constitutes the induction basis. To check the hypothesis for arbitrary k 3 , we assume that 
jF(/c' l5 k' 2 , k 3 ) is known for arbitrary k[ 2 and find a relation between !F(ki, k2, k 3 + 1) and 
T{k Xl k' 2 , k 3 ). Simple trigonometric transformations together with integrations by parts yield 

T{k u k 2 , k 3 + l) = T{k x + 1, k 2 , k 3 ) + T{k u k 2 + 1, k 3 ) - 2T{k x + 1, k 2 + 1, k 3 ) + (C3) 

2 



+ 



(An + 1)(A; 2 + 1) 



F®(k x + l,k 2 + l,k 3 ) 



where 

F ip \k u k 2 ,k 3 ) = 



(27T) 



2tt 



AS. 



sm 



2ki 



sm 



2A:_> 



sm 



2A:;j 



T {2 \Kk 2 M) = y 



(2A* - l)F(h - 1, fc 2 , h) - 2k x T{k x , k 2 , h) 
Inserting ( 1C5I) into ( 1C3f) . we get the relation 



F{k x , k 2 , k 3 + 1) = JF^ + 1, fc 2 , fc 3 ) + ( I : 

h + 1 



2fci + l\ 



-2 1 



fc 2 + 1 J 
F{k x + l,k 2 + l,k 3 



T{k x , k 2 + 1, k 3 ) 



,(C4) 
(C5) 

(C6) 



We substitute Eq. (1C2l) in the right-hand side of Eq. (1C6j) and derive the answer 



7T 



1/2 



71-1/2 7I-V 2 

r(A;i + A; 2 + fc 3 + 2) 



(C7) 



r(A; 1 + £; 2 + 1)T(A; 1 + k 3 + 2)T(k 2 + k 3 + 2) ' 
which satisfies Eq.( 1C2l) . Thus the induction over k 3 is verified and Eq.( 1C2l) is proven. 
Using flC2IIC4IIC5j) . one can also show that 



F ( ? n ~ 1 \kx, k 2 , k 3 ) = 0, n = l,2,... 



F {2 \k x ,k 2 ,k 3 ) = T{k x ,k 2 ,k 3 



k\k 2 k 3 



k\ + k 2 + k 3 



(C8) 
(C9) 



and 



F {4 \h,k 2 ,k 3 ) 



2(k x - l)(2k x - l)(2k t - 3)^(^1 - 2, fc 2 , A; 3 



-4(4fc? _ Qkl + Akx - l)T{k x - 1, k 2 , k 3 ) + 8klF(k h k 2 , k 3 



-^ 2 \k u k 2 ,k 3 )(l 



( fcl -i)(fc 2 -l)(fc 3 -l) 
&i + & 2 + & 3 — 1 



(CIO) 



APPENDIX D: AVERAGING Str fc [QiQ 2 ] AND Str fcl [QiQ 2 ]Str fc2 [QiQ 3 ]Str fc3 [Q 2 Q3] 
OVER THE PHASES AT i?i j2 , 3 = AND INTEGERS POWERS k, k 1<2 , 3 



over the phases 0i 2 : 

-Rl,2=0 



Let us average S± 2 = (Str[QiQ 2 ]) fc 

SiS 2 ( . 2 (<f>\ . n( * * \ • 2 f (P\ . cos (0) * * ^ 
— ^— I 4 sin I - J + 2(a A a A - a^a^) sm I - I H — a A a A a J? aij 1 



OK 
°12 



(Dl) 



where is integer; = 12 + A 12 ; 12 = <pi — 2 and indices 1 and 2 mark the variables 
of the matrices Q\ and Q 2 , respectively (see also notations in the Appendix We have 
to calculate the following integral of the periodic function: 

"1,2^2(01- 02 + A 12 ) = — / d0 12 «Sf 2 (0 12 ). (D2) 



(2tt) 2 77 " rM " Jivn " z ' 1Z/ 2tt 7 

The nilpotents from A 12 = ~ (A^ 2) + A^ 12) ), A^ = {t) R /a)% (vr/a)i ~ c.c, give no contri 
bution to the integral ( 1D2I) due to the periodicity of the integrand. 

Collecting the terms with the same powers of the Grassmann variables we find 

SiS 2 \ k f/„ . ,^\ 2k k/ . . .\2* 



5 * a= l^rJ H 2sin(0/2) i +2v 2sin(0/2) J («W -«£«*)+ (D3) 

2(fc-i) k ( , , A 2fc " 



fc 

+ 2 



2 sin (0/2) J --(2 sin (0/2)) 



a A a A a R a R 



We insert Eq. flD3p into the integral flD2p and integrate over the phase i2 using Eq. (ICl|) : 



27T / c C \ k 



's 1 s 2 y r(2fc-i 
2vr 7 -^-^^ v 4 / r(fc- i)r(fc) 



12^(012) = ( 1 w v~" x (D4) 



^ _ ^ (^2 + k(a* A a A - «#«??)) - fca^a^a^ajj | . 
Note that the equality ( 1D3I) can be rewritten in a more compact form: 

S f 2 = (S A f {(-n<*/2))V» - ^(si„ W2)) 2M T? 2 } ;T H S . 

(D5) 

Unlike Eq.( 1D2l) written for the case of 2 linked supermatrices, the averaged product 
<Sii <S 2 i (which includes 3 linked supermatrices) depends on the nilpotents coming from 
A® = A 12 + A 13 + A 23 and it can be expanded in the even powers of A^: 

ii,2,3 (0! - 2 + A 12 ) 5g (0i - 3 + A 13 ) S2I (02 - 03 + A 23 ) = (D6) 



(27T) 







j±p ff% d0' (0) Sll (0') S£ (0 - 0' + A( 3 )) 



* (A( g )) 2p 

h (2p)! 



(2vr) 2 



2tt 



d0 d0' (0) (00 $*S£ (0 - 00 



(D7) 



Here fci )2i 3 are integer. There are no odd powers of A® due to the equality (|(J8jl . For the 
purpose of the present paper, we need only the terms with p — 0,1, 2. The term oc (A®) 6 
yields zero after multiplication by the factor lZ p A q (see Eqs.( )39ll38i) ). Using the formula 
(ID5I) and the results of Appendix O, we find: 



(A( 3 )) 2p 



(2tt)- ,,„ 



2;t 



d0 d0' s£ (0) sg (00 dps* (0 - 00 
(a( 3 ))' 



p=2 



4! 



(D8) 



(A( 3 )) 



2p 



(2 



^ JJ% d0' Sg (0) Sg (00 ^Sg (0 - 00 



p=l 



(D9) 



(5 1 5 2 ) fcl (5i5 3 ) fc2 (5 2 5 3 ) fc3 -^P-(2^ 2 )(A; 1 , fc 2) fc 3 ) {hT 12 + A; 2 T 13 + fc 3 T 23 ) 2 - 



- [^ 2 )(A:i - 1, k 2 , k 3 )hrl 2 + J&\k u h - 1, fc 3 )A; 2 T 2 3 + J&\ki, k 2 , k 3 - l)k 3 T 2 23 ] ; 



where 



( A (3)) 2P 



(27T) 



2jt 



d0 d0' Sg (0) Sg (00 <9M| (0 - 00 



p=0 



(S 1 S 2 ) kl {S l S 3 ) k2 (S 2 S 3 ) k3 (PaH 1 + Part 2 + Part 3 ) ; 



Parti 



(2 k x k 2 k 3 T 1j2 T 1i3 T 253 [A;iTi >2 + k 2 T h3 + k 3 T 2 , 3 } + 



+(k 1 T lt2 k 2 T 1>3 ) 2 + (fciT^T^) 2 + (* 2 T 1>3 A;3T 2 • 



(D10) 



Part 2 = -~(r(k 1 -i,k 2 ,k 3 )k 1 rl 2 (k 2 r h3 + k 3 r 2!3 ) 2 + 

+F(k u k 2 - 1, k 3 ) k 2 Tj 3 (A; 1 T lj2 + A; 3 T 2i3 ) 2 + T{k x , k 2 , k 3 - l)k 3 T 2 2 3 (A; 1 T li2 + k 2 T x 



(Dll) 



3; , 



Part 3 = i (^(A* - 1, A; 2 - 1, fc-^T 2 2 £; 2 r 2 3 + (D12) 
+F{k x - 1, k 2 , k 3 - l)hTl 2 k 3 Tl 3 + F(k u k 2 -l,k 3 - l)k 2 Tl 3 k 3 T 2 2 ^ . 



APPENDIX E: RESULTS OF THE INTEGRATION 
OVER THE GRASSMANN VARIABLES 



A direct integration over the large number of the Grassmann variables (12 Grassmanns in 
the case of 3-matrix approximation) is technically trivial but very long and boring arithmetic 
procedure. We have used the "Grassmann" package of the "Maple" system to do this step 
of the calculations. Here we present the results of this procedure which are necessary for 
the calculations in the 3-matrix approximation 



4' 

d{V* R VRVlVA}p,q,m {V*RVR) p (vlVA) q (^T = ^ 



(El) 
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d{v* R VRVlVA} P , q ,m {v*rVr) AvIVa) a (A®) 2 x 



T 

p,q 

p,m 

•y-2 

9> m 



> = o 



J d{ri* R VRVAVA} P , q ,m (v r Vr) p {VaVa) q (A (3) ) 2 x 



p,q -L q,m 
J- p,m -L m,q 
q,p J- p.m 
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(E2) 



(E3) 



/ dW R VRVAVA}p,g,m {v*rVr) p (v*aVa) q X < 



•y-2 y2 

p,q q,m 

y2 y-2 

p,m m,q 

y2 y-2 

1iP P> m 



f o2 ' 



/ d {v R VRVAVA} P , q , m {vrVr) p (vaVa) ( 



x < 



^2 r y n y 

1 p,q 1 q,m 1 m,P 

nyn nyn2 nyn 

1 P,<? 1 g,m 1 m,p 

nyn nyn nyn 2 

1 p,q 1 g,m 1 m> p 



1 



(E4) 



(E5) 



see the definitions of the nilpotents in the previous Appendix. Note that all 3 indices p, q 
and m are different. 
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